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Abstract: We show that several types of global strings occur in colour superconducting 
quark matter due to the spontaneous violation of relevant U{1) symmetries. These include 
the baryon U(1)b, and approximate axial U(1)a symmetries as well as an approximate 
U(l)s arising from kaon condensation. We discuss some general properties of these strings 
and their interactions. In particular, we demonstrate that the U(1)a strings behave as 
superconducting strings. We draw some parallels between these strings and global cos- 
mological strings and discuss some possible implications of these strings to the physics in 
neutron star cores. 
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1. Introduction 

Domain walls and strings are common examples of topological defects which are present 
in various field theories [1, 2]. Domain walls are configurations of fields related to the 
nontrivial mapping ttq(M), while topologically stable strings are due to the mapping tti(M) 
where M is manifold of degenerate vacuum states. While it is generally believed that 
there are no domain walls or strings in the Standard Model due to the triviality of the 
corresponding mappings, such objects may exist in extensions to the Standard Model, or 
in phases where the symmetries of the Standard Model are broken. 

Topological defects that occur in extensions to the Standard Model may play an im- 
portant role in cosmology as described in the book [2], however, the focus of this paper will 
be strings and domain walls that exist in high density matter where the symmetries of the 
Standard Model are broken. As was recently demonstrated [3], in the regime of high baryon 
densities when the chemical potential fi is much larger than the QCD scale fi 3> Aqcd ; 
QCD supports domain walls. The existence and long life time of these domain walls is 
based on the following facts: 1) the instanton density is small at large chemical potential, 
suppressing the effect of the chiral anomaly and giving rise to the approximate U{\)a sym- 
metry; 2) this U(1)a symmetry is spontaneously broken, and 3) the decay constant of the 
pseudoscalar singlet boson (7/) is large and its mass small at large \i. 

In this paper we show that, along with these domain walls, there exist topological string 
configurations with interesting physical properties. Specifically, we discuss the three flavour 
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Nf = 3 colour-flavour locking phase (CFL), and the two flavour Nf = 2 superconducting 
phase (2SC). In the Nf = 3 CFL phase both the U(1)b and the approximate U(1)a 
symmetries of QCD are spontaneously broken. In addition, a third approximate U(l)s 
symmetry, related to the formation of a kaon condensate, is broken (see Section \2.2\ below). 
Thus, at least three types of strings are possible in the presence of a kaon condensate 
and two types of strings are possible if no kaon condensate forms. In the Nf = 2 2SC 
phase, the U(1)b symmetry is not broken and we only expect one type of string from the 
breaking of the approximate U (1)a symmetry. The important point is that the asymptotic 
freedom of QCD allows us to assert the existence of these strings in the high baryon density 
regime, and the properties of these strings can be determined by controllable weak-coupling 
calculations. 

Naively, these strings are neutral objects with respect to electromagnetic interactions, 
and therefore, one might think that they play no role in electromagnetic dynamics. How- 
ever, for the U(1)a strings, this naive expectation is incorrect due to the anomalous coupling 
of the singlet Goldstone field with a massless "photon" — actually a mixture of the bare 
electromagnetic field and one component of the colour gluon fields. We demonstrate that 
such strings behave as superconducting strings, and we estimate the maximum current 
that they can carry. Thus, QCD strings at high density may affect the electromagnetic 
properties of high density matter, influencing, for example, the magnetic fields in the cores 
of neutron stars. 

Regarding the baryonic U(1)b strings present in the CFL phase: it has been known 
for quite a while [4] that, although superfluid vortices (in liquid helium for example) and 
static global string solutions (related to the spontaneous breaking of a global symmetry, 
U(1)b in this case) are closely related, there are some differences. In particular, fluid 
vortices carries angular momentum while a static global strings have no angular momentum. 
What was demonstrated in [4] was that the two can be identified if the static strings are 
immersed in a uniform background with nonzero density. In high density QCD, such a 
background is naturally present, thus we demonstrate that the U(1)b global strings can be 
identified with the superfluid vortices likely to form in a rotating CFL phase with nonzero 
angular momentum as has been conjecture to exist in the cores of neutron stars (see for 
example [5, 6]). A similar mechanism leads to the formation of U(l)s strings. 

The mechanism that forms the U{V)a strings in the 2SC phase, which do not carry 
angular momentum, is less obvious to us and, in this paper, we simply assume that some 
mechanism exists to form these types of objects. In general, there will be complicated 
dynamical interactions between the strings, which may lead to the formation of more 
complicated stable objects like rings, springs and vortons. In this paper we limit ourselves 
to describing the properties of individual strings and their dominant pairwise interactions. 

This paper is organized as follows: In Section ^ we describe the string solutions and 
calculate the relevant parameters of the string in terms of the fundamental parameters 
of the theory. In Section || we discuss how the U(1)b strings present in the CFL phase 
are not static, but rather spin with a time dependent phase. These spinning strings can 
then be identified with the superfluid vortices. In Section || we discuss the anomalous 
electromagnetic properties of the U(1)a strings. We solve the relevant Maxwell equations 
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in the string background and demonstrate that the U{1)a strings behave as if they are 
superconducting. Finally, in the Conclusion ||, we discuss the interactions between strings 
and speculate on their effects on the physics of high density quark matter and neutron 
stars. 

2. Strings in High Density QCD. 

When one considers QCD at high density, the relevant excitations are due to quarks with 
momenta close to the Fermi surface. At high densities, these low-energy excitations have 
high momenta and, because of asymptotic freedom, one might hope to regain theoretical 
control in a weak coupling regime. Although the situation is not quite as simple as it 
appears, we do gain theoretical control. In particular, the colour anti-triplet 3 single gluon 
interaction is attractive, and, for low energies, this leads to a preferred Bardeen-Cooper- 
Schrieffer (BCS) pairing of quarks with opposite momenta that reduces the energy of the 
vacuum state. This is the phenomena of colour superconductivity. (For a nice review, 
see [7].) 

To be specific, we consider the simplest model where high-density U(1)a type strings 
appear: QCD with Nf = 2 massless quark flavours {u and d) and N c = 3 colours. This 
model is a rather good approximation to realistic quark matter at moderate densities. At 
higher densities, the approximation of Nf = 3 massless quarks becomes quit good. This will 
be discussed later in Section ^L^. The most important qualitative difference between the 
Nf = 2 and the Nf = 3 phase is the emergence of new spontaneously broken symmetries: 
the U(1)b in the three flavour case and the U(l)s symmetry if a kaon condensate forms 
in the CFL phase. As a consequence, if a kaon condensate forms in the CFL phase, then 
there emerge two new types of global strings related to the U(1)b and U(l)s symmetries. 
If a kaon condensate does not form, then only one new type of string related to the U(1)b 
symmetry emerges. In what follows we use the same normalization factors as the paper [3] 
on domain walls in dense QCD. 

2.1 Nf = 2 Superconducting Phase (2SC) 

We recall that the ground state at high baryon densities is a superconducting state [8, 9, 
10, 11, 12, 13], characterized by the condensation of diquark Cooper pairs. The super- 
conducting ground state spontaneously breaks the symmetry of the bare QCD Lagrangian 
through the non-zero diquark condensates ) which represent the Cooper pairs. Here 

we explicitly show the flavour (a, (3, etc.), colour (a, b, etc.) and spinor (i, j, etc.) indices. 
In this section we consider the Nf = 2 case where the strange quark is treated as heavy. 
In this 2SC phase, the diquark condensates have the form: 

Wa*f 6 >* = ^e^e abc X c , (2.1a) 
(^ l ^ b r = e^e^e abc Y c . (2.1b) 

The condensates X c and Y c are complex colour 3- vectors which are aligned along the same 
direction in the ground state. They spontaneously break the colour SU (3) c group down 
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to SU(2) C . The lengths of these vectors are equal \X\ = \Y\ and have been computed 
perturbatively [3, 14]: 

\X\ = \Y\=*£*. (2.2) 
2V2-JT g 

In perturbation theory, there is an approximate degeneracy of the ground state with respect 
to the relative U(1)a phase between X c and Y c which is a symmetry of the QCD Lagrangian 
at the classical level. A nonzero value (^4|) for the vacuum condensate implies that the 
U(1)a symmetry is spontaneously broken and, thus, the corresponding pseudo-Goldstone 
boson — the rj — enters into the theory. The U(1)b symmetry of the QCD Lagrangian also 
appears to be broken, but is in fact restored by a simultaneous SU(3) C rotation. Thus, 
only the U(1)a symmetry is spontaneously broken in the Nf = 2 2SC phase. 

Following [3], to describe the rf physics in an explicit way, we construct the gauge- 
invariant order parameter 

£ = Y ] X = Y*X C , (2.3) 
where the phase ifA is to be identified with dynamical r/ field, 

£ = |S|e-^ A = |£|e- i7? '//. (2.4) 
It is evident that £ carries a non-zero U(1)a charge: 

^ ^ e i75«/2^ ^ x ^ e -ia &nd Y ^ g ia ^ ^ 

£ -» e~ 2ia , (2.5b) 
=> ifA — ► <Pa + 2a. (2.5c) 

At low energies, the dynamics of the Goldstone mode <pa are described by an effective 
Lagrangian, which, to leading order in derivatives, must take the following form [3], 

C = f 2 [(do<p A ) 2 - u 2 {d.^ A ) 2 ] - V inst (ip A )- (2.6) 

For large chemical potentials fi 3> Aqcd> the leading perturbative values for the decay con- 
stant / and velocity u [15], and the instanton contribution describing the explicit anomalous 
breaking of the U{1)a symmetry [3] have been calculated: 

V inst ((pA) = -an 2 A 2 cos (p A , (2.7a) 

f2 = J*!_ u 2 = I (2.7b) 
1 8tt 2 ' 3" V ' 

In this formula A is the BCS gap, and a is a dimensionless parameter [3] 



29/3 

a ^[^±±\ (2.8) 



Aqcd^ 
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that goes to zero at large p. In what follows it will be important that the rf mass is 
asymptotically small, as can be seen from ( |2.6| ): 



(2.9) 



The effective Lagrangian ( [2.6D is justified for describing the light degree of freedom, 
but to describe global strings we must formulate an effective theory for fluctuations in the 
magnitude of the condensate From (|2.2| ) we have that = 9^ 4 A 2 /(87r 2 g 2 ). Thus, 
we introduce the dynamical field <£(x) of dimension one and expectation value \(&)\ = 
(p)=A: 



E(x) = A 



With this definition, the rj' dynamical field is merely the phase of a complex field 

= |$| exp(— if a) = /oexp(— if a) 



(2.10) 



(2.11) 



as in the Abelian Higgs model. 

The effective potential for this condensate has been calculated [16] for asymptotically 
large p. In terms of the field p{x) = \§(x)\, an approximate expression 1 for the potential [16] 
can be used which is a good description for p close to its vacuum expectation value p « 
(p) = A 
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p 2 A 

2~P 
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(2.12) 



Finally, in terms of a single complex field = pexp(— if a) the effective Lagrangian de- 
scribing the rf phase ipA and the absolute value for the condensate p can be represented in 
the following simple way: 



C 



eff 



\do$\ 2 



2 ia<Di 2 ' 



' A 2 

^((dop) 2 -u 2 (d lP ) 2 



p 2 A 



7T- 
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+ ap Apcos ip a 



p?A 

7T Z 



A 2 

1 



(2.13) 



+ ap Apcos (p a 



where the normalization factor f 2 /A 2 for the kinetic term has been chosen to repro- 
duce ( p]| ) where p ~ (p) which correctly describes dynamics of the light pseudo-Goldstone 
r/ meson. 

We should comment here that the potential presented in ( [2,12 ) was derived in [16] 
for very large p and when the field $ is close to its vacuum expectation value {p) ~ A. 
It deviates considerably from this form when p is far from its vacuum expectation value 



1 Note, the normalization for the condensate in [16] differs from our normalization for \X\ by a factor of 



16. 
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p ^ A. Besides that, there is no justification to keep only the lowest derivative term for 
the massive mode p in the expression (^]l]). Finally, there is an ambiguity in the definition 
of the dynamical field p describing heavy ~ A degrees of freedom: any smooth function of 
F{p) (for example exp(/o)) is appropriate for the description of the dynamics of the heavy 
degrees of freedom. This is a marked contrast with description of the light Goldstone fields 
where the physics does not depend on specific parameterization of the light fields. There 
are many other deficiencies of the Lagrangian describing massive field p which we 

shall not comment about. 

We are not pretending to have derived a Lagrangian describing a heavy (order A) 
degree of freedom. Rather, we want to demonstrate the qualitative features of the effective 
potential: that it has a Mexican hat shape (as it should) and that the internal relevant 
scales are of order A and not p. Indeed, the general scale ~ p 2 factors from the expres- 
sion (|2.1| ). Therefore, the description of the strings which follows serves mainly to illustrate 
the qualitative features of the strings. In particular, the details of the strings core are not 
well founded, however, we shall see that the behaviour of the strings far from the core is 
governed by the form of the potential where p ~ (p) = A where the effective potential is 
valid. This region gives a logarithmic contribution to the string tension which dominates 



the energy of the strings. Thus, despite the fact that the effective theory (2.1) breaks down 
for small p, the description of the strings far from the core is well justified. 

2.1.1 Global Strings 

First we consider the properties of an isolated global string which is symmetric about the 
z axis. The term Vi ns t((fA), which explicitly breaks the U(1)a symmetry, is small, and 
vanishes in the high density limit. Thus, our approximation of the global string is a good 
description out to lengths scales of order rriZ . 

With this simplification, we are looking for a static, classical field configuration <£(r, <p) = 
p(r, 4>) exp(—iipA(r, 4>)) which minimizes the energy density or string tension a: 

a = jj (n(p, <p A ) - H vac ) dxdy (2.14) 

where 7i is the energy density of the field configuration: 



1 - In 



and Wvac = T~l((p), {va}) = — p 2 A 2 /ir 2 is a trivial shift of the background vacuum energy 
introduced so that Ti — > far from the core of the string. 

To simplify the equations, we introduce the dimensionless field p and the dimensionless 
coordinates x^. 

p = (p) = 1, Xi = mxi = 2V3Axi. (2.15) 

The natural length scale is thus set by the parameter m = 2y3A: the mass of the excita- 
tions about the condensate p = {p)+5p in our model fl2.12|) . In terms of these dimensionless 
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parameters, the energy density becomes 

H = 24/VA 2 (i(0 4/ 5) 2 + \p 2 {d^ A f - p [1 - lnp] 

where all the derivatives are with respect to the dimensionless coordinates 5j. To minimize 
the string tension, we can drop the overall factor and then determine the equations of 
motion. To achieve the appropriate boundary conditions, ipA wm wind uniformly n times 
as a function of <ft 

<p A (r,<j>) = ¥a{4>) =n(p = ntan" 1 . (2.16) 
Converting to polar coordinates we have 

a = 47T/V J* (^ + ^- + V{p)^ fdf (2.17) 

where the dot signifies differentiation with respect to f, 

V(p) = l-p[l-lnp], (2.18) 

and we have introduced an outer limit R to the string's size to make the tension finite. 
The equations of motion follow from an application of a variational principle: 



P = pn 2 dV(p ) 
f f 2 dp 



with boundary conditions p(0) = and p(oo) = 1. The solution for the lowest energy 
string with winding n = 1 is presented in Figure |]. 

The string is governed by two relevant parameters: its core size and its tension. As 
can be seen in Figure p], the core size is of order r c ~ 1 while the tension of a global 
string diverges as log(-R) where R is an upper cutoff determined by the environment of the 
string. For large distances the second term in ( p. 17 ) dominates and we have a logarithmic 



divergence. We plot the cumulative energy as a function of the upper cutoff R in Figure |1[ 
Physically, the effective potential ( p. 12 ) is only valid for p ~ (p). Thus, the details of 



the string core must be interpreted with caution. The large distance behaviour and the 
logarithmic divergence in the string tension, however, are well justified. 

2.1.2 Domain Walls 

The formation of global strings discussed in the previous section neglected the instanton 
contribution Vi ns t((fiA) which explicitly breaks the U(1) A symmetry responsible for the 
global strings. Thus, the preceding analysis and picture is really only justified on distance 
scales small compared to those set by the anisotropy Vi nst (tfA), i.e. for r <C R ~ m~} ', 
which, in the high density limit, is much larger than the core size A as can be seen from ( |2.£| ) 
and (|2.9| ) . For distances larger than r 3> R ~ rr \^ > the appropriate description is no longer 
one of global strings, but one of QCD domain walls bounded by strings: the situation is 
similar to the so-called N = 1 axion model [17]. 
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Figure 1: Radial p(f) dependence of the n = 1 string in dimensionless units: the vertical scale is 
set by A and the horizontal scale is set by (2\/3A) _1 . 

Here we summarize a few results regarding QCD domain walls that will be relevant 
for our discussions later on. As classically stable objects, QCD domain walls were first 
analyzed in [18] for [i = in the large N c limit. Similar objects were later shown to be 
stable for N c = 3 at high densities [3]. We shall not repeat the analysis presented in [3] 
but quote some important results that we will make use of later. 

The thickness of the domain walls is set by the mass m v ' ( |2.9| ) of the excitations in the 
ifA field about the true vacuum (p^ = 0. In addition, the energy density per unit area of 
the domain walls (the wall tension) is 



Thus, on scales much smaller than R ~ m~, , the description of the global strings is 
valid, however, when one looks at scales larger than this, then one sees that these strings 
are really attached to domain walls. 

2.2 N f = 3 CFL Phase 

For very high densities, it can be a good approximation to neglect the strange quark mass. 
In this case there will be three flavours Nf = 3 in our model. Again, the ground state is a 
superconducting state, characterized by the condensation of diquark Cooper pairs [19, 20]. 
The superconducting ground state spontaneously breaks the symmetry of the bare QCD 



a = 8y2auf fj,A ~ ^/afJ 2 A. 



(2.20) 
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Outer string radius i? (in units of (2\/3A) _1 ) 

Figure 2: Cumulative contributions to the string tension of an n — 1 string from the core out to 
the cutoff radius R in dimensionless units: the vertical scale is set by AttJ 2 u 2 and the horizontal 



scale is set by (2V3A) -1 . The three curves correspond to the three terms in ( 2.17 ). The first kinetic 
term approaches 0.11 while the potential term approaches 0.25. Notice that the second term gives 
a logarithmic divergence in R. 



Lagrangian through the non-zero diquark condensates j^\Hj 6 ) which represent the Cooper 
pairs. (Again, we explicitly show the flavour (a, (3, etc.), colour (a, b, etc.) and spinor (i, 
j, etc.) indices). This condensate has the form 

<*ft*S>* ~ eVe^e abc X°, (2.21a) 
<*&*S>* ~ e^e^e ahc Y^. (2.21b) 

Now the condensates and are complex colour-flavour matrices. Following [21, 22], 
we introduce 



= ytjjf = (Y*)cXp. (2.22) 

The matrix S is a colour singlet and describes the meson octet as well as the ipA axial 
singlet. Unlike the 2S'C phase, the U(1)b symmetry is spontaneously broken giving rise to 
both axial and baryonic strings. 

To be precise, we consider the CFL phase and define the singlet phases if a and (fB 
describing the Goldstone bosons related to the spontaneously broken symmetries U{1)a 
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and U(1)b through the following structure 

<*&*S>* ~ e«e°^ W^ - *^ , (2.23a) 
<*&*S>* ~ e^'e^W*^"^- (2.23b) 

In terms of the </?a and ips fields, the corresponding strings can be described as was done 
above for 2SC phase: the only difference is that the U(1)b strings are not attached to 
domain walls. This is because the U(1)b symmetry is not explicitly broken: i.e. the 
equations describing the (fg string contain no terms analogous to those proportional to a 
in©- 

It was recently argued [23, 24, 25] that the so-called "symmetric CFL phase" — a pure 
form of the CFL phase where no other meson condensates are formed — is unlikely to occur 
in nature. Instead, it is likely that, in the CFL phase at high baryon density, a kaon 
condensate forms [23, 24, 25]. If this is the case, the kaon condensation breaks a new 
global U(l)s symmetry along with the U(1)a and U(1)b symmetries discussed above. 
From Goldstone's theorem one expects a new Goldstone boson tpg, which is the phase of 
the condensate, to appear in the spectrum. One can present the same arguments as before 
to deduce that a new type of the global string related to ips phase exists. We shall call this 
string the U(l)g string. The equation describing this new U(l)s string is the same as ( |2.1| ) 
with different values fs and as for the parameter / and coefficient a in the Lagrangian 
Q2.1| ). What is important is that, if the decay constant fs ~ [i is roughly the same as 
before and proportional to //, then the coefficient as is extremely small (but not zero) as 
was argued recently [26], 

a s ~ ^ 2 G F m u m s cos 6 C sin C , m^ s ~ a s A 2 . (2.24) 

In this formula Be is Cabibbo angle and Gp is the Fermi constant. Numerically, the scale 
corresponding to the pseudo-Goldstone boson mass m tps is of order m (ps ~ 50 KeV. 

3. Spinning Strings 

The main goal of this section is to demonstrate the connection between the idea of a global 
U(1)b string and the superfluid U(1)b rotational vortices in the CFL phase that carry 
angular momentum. Since they carry angular momentum, they may be copiously formed 
in a rotating neutron stars provided that the core of the star is in a colour superconducting 
phase. 

It is generally accepted that the global strings and non-relativistic superfluid vortices 
are closely related; however, the behaviour of the two are quite different. For example, 
near a straight superfluid vortex at rest there is a velocity field moving circularly around 
it, carrying momentum and kinetic energy. But a simple global string (for example defined 
by Lagrangian (|2.6|) ) is a time independent solution of the equations of motion. The 
momentum density away from the string core is ~ f 2 doipAdi<PA which is zero in the rest 
frame of a static global string. The precise relation between global strings which we discuss 
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in this paper and superfluid vortices was given in [4]. There it was demonstrated that, if 
a global string is immersed in a uniform, Lorentz-non-invariant background, then the 
static solution (with tps to be identified with the azimuthal angle) is replaced by a time- 
dependent solution (fB — » fB + wt where the coefficient u is determined by the density of 
the background and the magnitude of the condensate. 

In the quark superconducting phase, such a background is naturally present and the 
constant to is uniquely determined by the chemical potential \i. The relevant Goldstone 
phases if a and (fB are as defined in ( [2.23|) . The effective Lagrangian in the CFL phase for 



the axial Goldstone field <pa can be derived in the same manner as before (2.6) 



C A = flKdoiPA? - u\{dwA?] - V inst ( VA ) ■ (3.1) 

For large chemical potentials, \x 3> Aqcd, the leading perturbative values for Ja and ua 
have been calculated [21, 22, 27, 15]: 

Mnst(^) = -a'/i 2 A 2 cos(^) (3.2a) 
f\ = %A = \. (3.2b) 

The instanton contribution that explicitly violates the U(1)a symmetry has also been 
calculated [28] 

which again vanishes for large \i. 

To leading order in perturbation theory, the parameters fs and ub for the Goldstone 
boson associated with the spontaneous breaking of U(1)b symmetry are identical to the 
ones presented in ( p. 2b] ) [21, 22, 15]. However, the structure of the effective Lagrangian 
for the baryon Goldstone field is a little bit different due to the explicit presence of the 
chemical potential [i in the original QCD Lagrangian. 

In order to restore the dependence on fj,, one can use the following trick which, in the 
present context, was originally suggested in [29, 30] and consequently has been used in a 
number of papers [21, 22, 15, 31, 32, 33]. The idea is to make use of the fact that /i enters 
the QCD Lagrangian in the same way as the zeroth component of a gauge potential: 

= f (yfdu - m) * + /uf 1 "^, (3.4a) 
= * (iYd v + ^7° - m) . (3.4b) 

Therefore, one can formally promote the global U(1)b baryon symmetry to a local one 
by introducing a gauge field B u = (Bq,0) coupled to the baryon current with coupling 
constant \x: 

£qcd = * {i-fVv - m) * (3.5) 

where V u = d u + ifiB u . 
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Indeed, under the U(1)b rotations we have: 



* -► e ia #, (3.6a) 
<PB ^ + 2a, (3.6b) 
B v ^ B v — fi~ 1 d u a, (3.6c) 

which leave the microscopical QCD Lagrangian unchanged. An effective description must 
respect this symmetry, and therefore, in the effective Lagrangian description of the baryon 



Goldstone mode ips ( 2.23 ), one must replace the derivative d u <p a in (|3.1[) by the covariant 
derivative D v y>B = d u (pB + 2fJ-B u . In matter with uniform density, we fix Bq = 1 and 
Bi = so that 



and thus 



D ip B = {d ipB + 2/x) Di Lp B = di<p B (3.7) 



C B = f 2 B {{D^ B ) 2 -u 2 B {D^ B ) 2 ). (3.8) 



From this equation we see that we can restore the original form (without covariant deriva- 
tives) if the baryonic Goldstone mode (fs — the phase of the condensate ( |2.23| ) — receives a 
time dependence in the ground state 

<p B ~ e i2 ^. (3.9) 

Note, in the previous literature concerning the CFL phase [15, 21, 22], this time dependence 
was irrelevant. However, for the discussions of strings as emphasized in [4] this dependence 
is essential. In particular, following [4], to analyze QCD strings and their interactions 
in a nontrivial environment, one can introduce, an equivalent description in terms of a 
two-index antisymmetric tensor field B^ where 

fd^B ~ e^ uXa d u B x °. (3.10) 

Using this formalism, one can calculate the Lorentz force between strings, similar to the 
Magnus force in the non-relativistic limit. One can also introduce a quantity similar to the 
non-relativistic vorticity and demonstrate that it is quantized (in our notations) in units 
of /U . We shall not discuss these interesting topics in this paper: once the appropriate 
correspondence is made, the techniques of [4] can be applied. The important remark that 
we make is that, if the CFL phase with the global U(1)b symmetry is realized in the interior 
of a rotating neutron star where there is a non-zero chemical potential fi, then the global 
strings that form will be spinning and will carry angular momentum. Thus, drawing upon 
the analogy with liquid helium, we expect that, if the CFL phase rotates, then spinning 
global U(1)b strings will form. 

One consequence of this connection between angular momentum and the global string 
is that, if global strings are formed through the transfer of angular momentum, then there 
will be a correlation between the direction of the angular momentum and the sign of the 
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winding of the string. To see this, note that the angular momentum of the field ips is 
proportional to: 



M{j ~ J Xjd ipBdiif B - XidoipsdjipB- ( 3 -ll) 

The terms doips ~ fJ- pickup the sign of the chemical potential and the terms ditpB ~ n 
pickup the sign of the winding or topological "charge" n ( p. 16 ) of the string. In the core 
of a neutron star, for example, the angular momentum has a definite sign (as set by the 
rotation of the star) and the chemical potential has a definite sign (the core is composed of 
baryons, not antibaryons). Thus, the sign of the topological charge is correlated with the 
sign of the angular momentum. If the formation of these strings is related to the rotation 
of the bulk phase, then there will be an excess of one type of string (either positive or 
negative winding). 

Exactly the same arguments can be presented for the U(l)s strings mentioned at the 



end of the Section 2.2 which originate from kaon condensation. The only modification 
that must be made to Equation ( |3.8| ) is to replace the chemical potential \l by an induced 
chemical potential fj, e g ~ m1/2pp, see [23, 24, 25]. The argument based on ( 3.11 ) about 



the correlation between the sign of the topological charge of the string with the sign of the 
angular momentum, also holds. 

The situation with U(1)a strings is expected to be quite different: conservation of P 
parity implies that the number of strings and anti-strings must be the same. This is due to 
the fact that U(1)b strings do not transform into anti-strings under the exchange R ^ L 
while U(1)a strings do, as can be seen 2 from the definition ( 2,23j ). 



4. Electromagnetic Properties 

Up until this point, we have discussed the existence of global U(l) strings and their cor- 
respondence with superfluid vortices. However, since they only involve excitations close to 
the Fermi surface, there will not be enough to affect the thermodynamics of the supercon- 
ducting phases. In addition, the fields ifA and (ps are neutral, so one might naively expect 
that they have little effect on electromagnetic physics either. It turns out, however, that 
the axial strings <pa have non-trivial electromagnetic properties. 

4.1 Anomalous Electromagnetism of the rf 

In this section we are mainly concern with electromagnetic interaction of the neutral r/ me- 
son. The simplest way to derive the corresponding low-energy effective Lagrangian, which 
includes the massless electromagnetic field and the light rf field ( |2.6D , is to follow the 
standard procedure and consider the transformation properties of the path integral under 



the U(1)a chiral transformation (2.5). As is known, the measure is not invariant under 



these transformations due to the chiral anomaly: it receives an additional contribution 



2 We thank Misha Stephanov for presenting this argument to us. 
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5C = (a/2)d v Jp. The expression for the anomaly d v is well known and takes the form: 

2 2 N f 

dUJ » = \h N ^ da + h NcFp ^ Q2 ^ (41) 

where we have included the electromagnetic fields along with the gluon fields Cr" and 
their duals: 

Ffu/ = —£[j,u\<jF > G = —£^,uXaG . (4-2) 

One should note that the expression for the anomaly is an operator relation which is 
valid for any finite chemical potential /i. Indeed, the anomaly arises from the ultraviolet 
properties of the theory and is not sensitive to the finite chemical potential as long as the 
regulator fields are heavier than \i. 

The second step is an identification of the parameters from the Lagrangian 5L with 
the physical fields of the theory. For example, in QCD with \x = the electromagnetic 
field F^ v (^2) from the original theory is the observable physical field (in contrast with 
superconducting phases when /i ^ 0, see below). The parameter a of U{\)\ chiral trans- 
formation is identified with the physical rf field which is nothing but the singlet phase of 
the chiral condensate, 

{* R * L ) ~ e~ ia ~ exp (-i-^=\ . (4.3) 

Therefore, the anomalous effective Lagrangian of rf coupled to photons takes a familiar 
form: 

2 N f 



8Cw~~ = — ttN c FF V Q 2 f '-=. (4.4) 



In a similar manner, one can derive the well-known effective Lagrangian describing the 
famous 7T° — ► 77 decay. In this case, one should consider the expressions for the third 
component of the axial current, 

e 2 

^(«7m76« " d W d) = g^FF, (4.5) 

and identify the corresponding transformation parameter 03 with Goldstone mode tt° such 
that the effective Lagrangian takes the familiar form 

= ^FF^. (4.6) 

Now, we want to derive a similar expression for the effective Lagrangian describing 
interaction of the light rf field with electromagnetism in CFL and 2SC phases 3 As we 

3 Anomalous electromagnetic interaction for the SU(3) Goldstone modes in the CFL phase was discussed 
previously using a different approach [34]. As far as we know, the anomalous electromagnetic interaction 
with the singlet rf phase has not been previously discussed in the literature. 
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mentioned above, the operator expression for the anomaly (^4]) remains unchanged: the 
singlet phase a defined in ( |2.5| ) is identified with the physical Goldstone mode r/, which is 
now the phase of the diquark condensate instead of the chiral condensate 

S~e-^|S|, a = fA = ^L. (4.7) 

'' 2 2f A V ' 

This is not the end of the story however, because, as it has been known since [35], in dense 
QCD matter, the electromagnetic field strength F^ v and the electric charge e are not the 
appropriate physical quantities. Rather, a combination of the electromagnetic field 
with 8 th component of the gluon field A^ acts as a physical massless photon A field [35] : 



A — 


cos 9 + A 8 sin 9 


(4.8a) 


A 8 - 


-Afj, sin 9 + A* cos 9 


(4.8b) 


cos 9 = 


9 


(4.8c) 


\J g 2 + 4t] 2 e 2 


smO = 


2erj 


(4.8d) 


\J g 2 + 4t] 2 e 2 


e = 


eg 


(4.8e) 


\J g 2 + 4r] 2 e 2 



where the A is the physical "photon" field and e is the physical charge. The values for 
parameter r\ entering the expressions fl4.8|) are given by [35], 



noFL = ±, ™c = -±=. (4.9) 

We note that our expression for the angle 9 in terms of e and g is obtained from the 
one presented in [35] by changing g — * g/2. The difference is due to the non-standard 
definition of the strong coupling constant g in [35] (the absence of the factor 1/2 in front of 
g in the covariant derivative in [35] ) . We use the standard definition for the strong coupling 
constant such that the chiral anomaly is given by Q4.1| ). Using our normalization for the 
rf field and expressing the anomaly in terms of the physical electromagnetic field we 
arrive to the following effective Lagrangian describing the interaction of the rf with the 
electromagnetic fields in the CFL phase: 

(CFL) = _e^ ^ , m 

0L -^TT 4vr 2 2 ' [ 1 

It is interesting to note that half of this result is due to the original electromagnetic inter- 
action, ~ FF while the other half is due to the gluonic part of the anomaly ( |4.1[ ). One can 
repeat the same calculations for the 2SC phase with the result similar to ( 4.1Cp : 

2 

w-(2SC) e TjrfA (A 11) 

However, for the 2SC phase, only one sixth of this contribution is due to the original 
gluon term in the anomaly. One should also note that the same procedure determines the 
anomalous coupling constants ~ g^TQ describing the decay of a heavy gluon to an rj and 
a photon — the decay similar to p — > n'y in usual QCD. We shall not discuss further this 
physics involving the heavy particle (mg ~ A) in the present paper. 
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4.2 Superconducting strings 



In this section we wish to analyze the system of light particles rj (the phase if a) an d photon 
Tfxv Combining ( |3.1| ) and ( 4.1C| ) we have the following effective Lagrangian: 



CeB = f iiPWA? ~ u\d^ A f) - Mnst(^) " + ^ (if) ^ ■ 

(4.12) 

We approach this problem for large /x. In this case, as a first approximation, we can omit 
the instanton contribution VJnst as we did for the discussion of a single string. This allows 
us to treat the problem of string-electromagnetism exactly at large distances. If we turn 
on the small instanton term and explicitly break the U(1)a symmetry, then domain walls 
are also allowed and the problem becomes much more complicated. We speculate on a 
possible effects of these physically relevant domain walls in the conclusion. 



The effective theory (4.12) is almost identical to that studied in [36]. We begin our 
analysis by considering a single spinning global axial string lying along the z-axis and use 
cylindrical coordinates (r,(f>): 

$(f, r, (f)) = p(r)e^ ( *'* ) . (4.13) 

The radial solution p(r) was presented in Figure [l], but we will only be concerned with 
distances large with respect to the core size of the string r > r c ~ A. In this regime, the 
effective theory ( f4.12| ) becomes valid and we have the following equations of motion for the 
"electromagnetic" field: 

d ^ v = h*^^ = ( 4 - 14 ) 

where a = z 2 /4ir is the modified fine structure constant and the background field if a has 
the solution 

ip A (t, <p) = VAt + 4> where 4> = tan -1 — . (4-15) 

x 

We assume that the rate at which the string is spinning is small: 

VA = ^<: A" 1 . (4.16) 

The origin for this assumption is that, while ips strings have very high frequencies of 
order of fi, ( |3.9| ), the ipA strings do not spin on their own. They can, however, receive 
some angular momentum through the interaction with the surrounding medium an/or the 
(fB strings. The interaction of the Goldstone fields with other particles is expected to be 
suppressed by some power of fj, . Therefore, we expect ip A to be small with respect to all 
relevant scales. 
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The anomalous Maxwell's equations ( 4.12j ) in the presence of a time dependent global 



string background have the following form far from the core of the string (r ^> A 1 ) 

d£ r (ldB z dBA 2a 2a<p A 

— + ~ » = S * B r, (4.17a 

at \r oq> oz J irr it 

_§ + (i^)_I^) = ^ r _^ Bi , (4 . 17b) 

at \r Or r o<p J irr tt 
88$ dB z dB r 2a(p A 

-dt ~ -dr~ + -dz~ = (417C) 
l^) + ^ + l^ = 2a 

r or oz r o<p irr 

In the limit of (pA — * 0, these results reduce to those presented in [36] with the replacement 
2a — > a which is due to the extra factor 2 in our expression ( |4.10[ ) for the anomaly as 
compared with the axion model considered in [36]. 

We were not able to find a complete solution of these linear equations over all of space, 
however, a z-independent solution for radii larger than the core size but not too large 
A -1 ~ r c <C r <C (p^ 1 can be easily found: 

£ r = C+r - l+2a ^ + C-r- 1 - 20 ^, (4.18a) 
B^ = C + r- 1+2a l* - C^ 1 - 2 ^, (4.18b) 
B z = p A {C + r 2a ^ + C^r- 2a ^), (4.18c) 

with the other field components vanishing. Coefficients C± need to be determined by 
matching the solution ( |4.18| ) with the behaviour of fields in the core region r ~ A where our 
effective Lagrangian approach breaks down. In what follows, we estimate these coefficients 
using dimensional arguments. 

The solution ( 4.18] ) is reduced to the corresponding expression [36] in the limit (pA = 0. 



In this case, the solution becomes valid for arbitrarily large r. Our approximate solu- 
tion ( [4.18 ) is limited to the range A -1 <C r <C ¥>7 ■ At shorter ranges we must further 



develop our microscopic model. 

The most important feature of the solution ( [4.18 ) is that, for a given C + ^ or C_ 7^ 0, 



the solution has the Lorentz property Bj, = ±£ r — as if the string carries a light-like current 
four vector with current density along the z-axis j« = (j,0, 0, rbj) with the whole system 
spinning with small frequency if a *C 1. The presence of a small B z field in ( |4.18c| ) which 
is due to the spinning of the static solution can be easily understood by boosting to the 
local frame moving with the speed = —•pat at the point r. 

We should note here, that the presence of a current with the property ) tli = (j, 0, 0, ±j) 
in the system defined by (4.12) was established by analyzing the only large distance physics 



without formulating a microscopic model of the core. A similar result was obtained in [36]. 
There, however, the effective theory was taken to be exact, allowing for a detailed analysis of 
the core physics. In that case, it can be seen that such a current results from a charged Dirac 
fermionic zero mode present in the string background (see [37, 38, 39, 40] for more details). 
Unfortunately, in our case, the effective theory ( 4.12|) breaks down at short distance scales: 
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we only have theoretical control over the large distance physics, however we believe that 
the microscopical explanation of the property j M = (j,0, 0,±j) in our case is very similar 
to the microscopical explanation given in [36] : namely, that it is due to zero modes of the 
charged fermion field which travels inside the core of the string. 

The presence of zero modes (localized in the (xy) direction) in the string background 
is a very general property of such a background and is a trivial consequence of an index 
theorem, [37, 38, 39, 40]. Although this topic is beyond the scope of the present paper, we 
would like to note that, in the colour superconducting CFL or 2SC phases, the interaction 
between the gapped fermions close to the Fermi surface and the diquark condensate has a 
more complicated algebraic structure than the simple Yukawa coupling of a single fermion 
considered previously [38, 39, 36]. Nevertheless, one can explicitly demonstrate the presence 
of the fermionic zero modes which would have been the regular gapped excitations in the 
absence of the string background. Therefore, we believe that the microscopical explanation 
of the result j M = (j,0, 0, =tj) in our case is analogous to the explanation given in [36]. 
However, unlike the case in [36], we expect that both the coefficients C± are nonzero 
due to the lack of Lorentz invariance in our system. (The dispersion laws for in-medium 
fermions do not have a standard Lorentz- invariant form.) 

One of the consequences of localized charged fermionic zero modes in a string back- 
ground is superconductivity. Indeed, as was demonstrated in [38], the problem can be 
reduced to a two dimensional effective theory which describes massless charged fermions 
(our original zero modes in four dimensions), and photons (our physical "photon" fields 
An) which are coupled to the massless fermions through the physical charge e. It was also 
demonstrated in [38] that such a system describes a superconducting string in the sense 
that, if external electric field is applied along the z direction, it results in persistent current 
along the string. We shall not repeat all these well known arguments, which are based 
on the bosonization technique of the (effectively two-dimensional) localized zero modes. 
Instead, we refer the reader to the original paper [38]. 

Thus, we can only present a dimensional estimate at this time: For a given orientation 
of the string solution ipA (r, 4>) = +4> we expect that both C± of Equation ( 4.18; ) are of the 
same order and can be estimated dimensionally to be 

C_ ~ eA(A)~ 2a ^, C+ - eA(A) +2a /\ (4.19) 

This result behaves in the same way as if fermions trapped in transverse zero modes travel 
in opposite directions (without cancellation, C + / C_) at the speed of light. We should 
remark that, due to the linearity of Maxwell's equations, a linear superposition of two 
sources c+(j, 0, 0, +)) + c_Q, 0, 0, — j) produces a desirable superposition of the correspond- 



ing solutions (4.18). Such electromagnetic behaviour of the string is reminiscent of the 
behaviour of free superconducting strings [38, 41]. 

The question of how the external electric field may be generated is a different question 
which is not addressed in the present paper. However, we should remark here that, in the 
core of a neutron star, the current could be generated by the motion of the strings through 
external magnetic fields already present in neutron star [38]. 
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To conclude this section, we estimate the maximum current which can be achieved in 
a colour superconductor. Under an external electric field applied parallel to the axis of the 
string, the fermion current will grow. Eventually, however, it will saturate: If the mass 
of the fermions away from the string is A, then once the Fermi momentum of the charge 
carriers rises above pp > A, it will be energetically favourable for the fermion to leave 
the string. Therefore, one expects that the maximum current supported by a single string 
would be of the order 

eA 

jmax ~ — (4.20) 

which cannot exceed 2 • 10 3 A. 

We do not expect, however, that a single, infinitely long string is realized in nature. 
Rather, we expect that the strings will take the shape of a ring or organize into more com- 
plicated objects. We presently cannot say more about the results of the rather complicated 
dynamics of strings in dense quark matter in the presence of strong external magnetic 
fields. However, we refer to some of the results discussed in [38] and the textbook [2]. As a 
first step to understanding the dynamics of this complicated system, one has to understand 
the forces which act between the strings. We discuss this subject in the Conclusion. 



5. Conclusion 

In this paper, we have considered the existence of classically stable global string configu- 
rations in high density QCD. We have shown that, in the Nf = 3 colour flavour locking 
(CFL) phase, baryonic strings resulting from the U(1)b symmetry of the condensate can 
carry angular momentum. A similar situation arises for U(l)s strings if a kaon condensate 
forms. In addition, we have shown that the axial strings resulting from the approximate 
U(1)a symmetry may exist in both the Nf = 3 CFL phase and the Nf = 2 (2SC) phase. 

Even though the fields that form these global strings are neutral, we show that the 
U(1)a strings have non-trivial electromagnetic properties due to the anomalous coupling 
between the rj' field and a massless "photon" that is a mixture of the eighth gluon and 
photon fields. Due to this anomalous interaction, axial strings behave as superconducting 
string. 

As we have already mentioned, however, in general, strings will not be formed as 
infinitely long strings in isolation. The details of the interaction between strings is quite 
varied and complicated: stable structures like helices, rings and vortons may form, or 
the strings may form networks or tangle. Besides that, U(1)a strings as well as U(l)s 
strings will be bounded by domain walls which themselves can decay into string-antistring 
configurations. In order to understand the possible dynamics of the system one must know 
different forces which determine the dynamics. 

In this Conclusion we shall only comment on a few simple possibilities here and refer 
the reader to the textbook [2] for a discussion about general string dynamics. We shall 
consider the following interactions between strings: 

• Inter- Vortex Force: Fyy. 
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• Domain Wall Force: Fuyy. 

• Anomalous Electromagnetic Force: Fem- 

We shall briefly summarize these forces and provide qualitative estimates in terms of the 
parameters of the QCD phase. In the following we assume that two infinitely long global 
strings lie parallel to the z-axis with separation d. We label the strings by their topological 
"charge" n which corresponds to their winding number ( |2.16| ) with respect to the positive 
z axis. Strings of the same charge have the same orientation: String of opposite charge can 
annihilate. All forces are presented as a force per unit length of the string. 

Inter- Vortex Force The most important and well-known interaction is due to the ex- 
change of massless Goldstone particles. The force is estimated by considering the energy 
of a configuration of two strings. This approximation is valid for distances d much larger 
than the core size d » A -1 [42, 43] 

Fw ~ T — ^— • (5.1) 

This force is repulsive for strings with the same charge, and attractive for strings of opposite 
charge. Physically, strings of opposite charge tend to attract (energetically they try to 
annihilate and restore the vacuum everywhere) while strings of the same charge repel: 
strings of higher charge n > 1 will split into several strings of single winding which will 
tend to move away from each other. 

For baryonic U(1)b strings this result extends to arbitrarily large distances d because 
the Goldstone bosons are massless. For U(l)s strings, this description is only valid for 
distances d < m~ l ~ (50 KeV) -1 because the Goldstone bosons are in fact massive due 
to weak interaction [26]. A similar situation occurs for the axial U(1)a strings where this 
description is only valid for distances d < m~} because the Goldstone bosons is massive 
due to instanton effects ([O^). Beyond this range, this inter- vortex force falls away and the 
domain wall force starts to dominate. 

Domain Wall Force On distance scales d > m^ 1 , where m = m„i for U(1)a strings and 



m — m,p s ~ ^/asA for U(l)s strings, ( |2.24 ), one will see that the U(1)a and U(l)s strings 



are really embedded in domain walls of thickness m . The scale for this interaction is set 
by the corresponding domain wall tension a calculated for if a field in [3] and for ip$ field 
in [26]: 

Fdw ~ <t ~ f 2 m. (5.2) 

This force can be attractive for opposite charge strings if a domain wall connects them. 
Domain walls are not formed between strings of the same charge if we neglect the process 
of nucleation which is equivalent to the creation of the string-antistring pairs. 

Anomalous Electromagnetic Force For distances between A -1 <C d <C ip^ 1 the axial 
string behave as current carrying wires. Through the exchange of "photons" , these strings 
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thus interact with force 



jj a A 2 

Fem ~ 3 ~ ± y (5.3) 

a ir a 



This force is attractive for strings of the same charge and repulsive for strings of opposite 
charge. Physically, strings carrying current in the same direction want to bunch together 
into a "wire". In comparison with the other two forces, however, the anomalous electro- 
magnetic force is suppressed by a factors of a and A//i and can be neglected in comparison 
with these other forces. 

Besides that, when a U(1)a string moves at velocity v across an external magnetic 
field, the effective superconducting current should build at a very high rate, (see [38]). In 
such a situation, an electric field is also generated leading to the creation of particles, a 
process which we do not discuss here. 

One more effect which we would like to mention here and which deserves for a further 



study is the binding strings by the domain wall force (5.2). In a situation similar to the 
N = 1 axion model [17] such a configuration decays very quickly. In our case when U{\)a 
string is superconducting and can form a ring with a dragged magnetic flux crossing the 
ring, the situation may not be so obvious. One should better understand the system of 
superconducting string with attached domain wall in the dense background in the presence 
of a nonzero magnetic field, before one can make any conclusions regarding the system. 
We suspect that due to the interactions which were not present in the axion case, some 
stable configurations (like vortons) will be possible. 

The most likely place to encounter bulk high-density superconducting quark matter is 
in the core of neutron stars. Typically, neutron stars are rapidly spinning, and we suspect 
that during the formation of the core, this rotation would impart some angular momentum 
to the superconducting core. As we showed in Section [3L spinning global U(1)b and U(l)s 
strings carry angular momentum that is correlated with the string's charge. 

It remains to be seen whether these strings can have any impact on observable effects 
(they might, for example, affect glitches: sudden increases in the rotation frequency uj of 
neutron stars by as much as Alo/lo ~ 10~ 6 , or the magnetic field structure and evolution). 

There will likely be many other effects due to the interactions and electromagnetic 
properties of global U(l) strings that will play a role in the physics of superconducting 
phases of high density quark matter. These effects might be closely related to analogous 
effects studied for cosmic strings (we refer the reader to the textbook [2]), and studied in 
condensed matter physics. We hope that the present paper will initiate some activity in 
this direction. 
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